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I.  Background  and  Objectives 

Electrons  produced  in  the  bow  shock  of  missiles  can  play  an  important 
role  in  determining  the  chemical  and  physical  properties  of  hypersonic  flow- 
fields.  Due  to  non-equilibrium  overshoots,  electron  temperatures  can  rise 
above  10,000  K  in  such  flowfields.  Furthermore,  as  flight  velocities  increase, 
electron  number  densities  and  collisions  become  more  significant.  However, 
electron  collision  data  needed  in  simulations  of  these  flowfields  are  very  frag¬ 
mentary  and  often  not  available.  This  is  particularly  true  for  near-threshold 
excitation  processes  where  the  incident  electrons  barely  have  enough  energy 
to  drive  the  excitation  and  for  collisions  with  metastable  states  and  with 
fragments  such  as  OH. 

Our  strategy  for  obtaining  these  electron-collision  cross  section  data,  par¬ 
ticularly  those  that  are  experimentally  inaccessible,  has  been  to  calculate 
them.  Calculations  of  electron-collision  cross  sections  at  the  low  energies  of 
interest  in  flowfields,  are  also  difficult.  In  fact,  although  the  physical  principles 
in  low-energy  electron-molecule  collisions  are  well  understood,  and  although 
several  methods  have  been  developed  for  numerical  studies  of  these  collisions, 
progress  in  their  application  has  been  very  limited.  What  is  required  is  not 
just  a  method  and  algorithm  by  which  these  cross  sections  can  be  calculated 
in  principle,  but  an  implementation  that  makes  relevant  problems  feasible. 
The  thrust  of  our  effort  has  hence  been  to  develop  innovative  scalable  imple¬ 
mentations  of  our  formulation  of  electron- molecule  collisions1,2  with  which  we 
could  harness  the  computational  power  of  larger  parallel  computers  to  obtain 
electron-collision  data  needed  in  modeling  flows. 

II.  Summary  of  Progress  and  Accomplishments 

At  the  low  energies  of  interest  in  these  flowfields,  an  accurate 
quantum-mechanical  treatment  of  the  full  3(N+l)-dimensional  space  of  the 
Schrodinger’s  equation  that  governs  the  collision  of  an  electron  with  an  N- 
electron  molecule  is  necessary;  low-order  approximations,  such  as  the  Born 
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approximation,  are  not  applicable.  Furthermore,  because  direct  numerical 
integration  of  Schrodinger’s  equation  in  this  3(N+l)-dimensional  space  is  not 
practical,  most  recent  studies  rely  on  the  use  of  variational  approximations  to 
obtain  the  scattering  amplitudes  for  these  collisions.  Our  studies  are  based 
on  the  Schwinger  multichannel  (SMC)  method,1  an  extension  of  Schwinger’s 
variational  principle.  This  SMC  variational  principle  was  specifically  devel¬ 
oped  for  applications  to  low-energy  electron-molecule  collisions  and  possesses 
compelling  computational  advantages.  However,  these  come  at  the  cost  of 
evaluating  a  Green’s  operator  whose  matrix  representation  is  obtainable  only 
by  numerical  quadrature  of  a  class  of  integrals.2  Indeed,  evaluating  and  trans¬ 
forming  the  data  for  this  quadrature  is  the  principal  computational  task  in 
applications.  Not  surprisingly,  this  task  had  been  a  bottleneck  and  hindrance 
to  wider  applications  of  the  SMC  on  sequential  computers  such  as  the  Cray 
Y-MP. 

We  have  successfully  developed  efficient,  scalable  algorithms  for  carrying 
out  the  main  computational  tasks  in  these  calculations  on  large  distributed- 
memory  parallel  computers  consisting  of  several  hundreds  of  interconnected 
microprocessors.  Our  strategy  for  “parallelization”  of  these  two  main  tasks 
—  namely,  evaluation  of  a  large  number  (typically  1012  to  1015)  of  elementary 
integrals  and  their  subsequent  transformation  to  obtain  the  physical  matrix 
elements  —  is  as  follows.  The  first  of  these  tasks  is  perfectly  parallel,  mean¬ 
ing  that  on  a  multiprocessor  machine  different  elementary  integrals  can  be 
evaluated  simultaneously  on  different  processors  with  effectively  zero  paral¬ 
lel  overhead.  Although  the  second  step  —  transformation  of  integrals  stored 
in  the  microprocessors  —  involves  communication  and  synchronization  over¬ 
heads,  it  can  be  formulated  in  terms  of  multiplication  of  large,  distributed 
matrices,  a  highly  efficient  procedure  on  parallel  computers.  With  these  scal¬ 
able  algorithms  we  routinely  achieve  aggregate  performance  in  the  tens  of 
gigaflops  range.  This  capability,  which  to  our  knowledge  is  unmatched,  has 
allowed  us  to  carry  out  sophisticated  first-principles  calculations  of  electron 
collisions  with  species  such  as  N2,  CO,  NO,  H2O,  OH,  and  CO2. 


2 


Highlights  of  these  application  include 

•  In  a  joint  effort  with  the  experimental  group  of  Professor  H.  Ehrhardt 
of  Kaiserslautern  University  in  Germany,  we  have  completed  a  comprehen¬ 
sive  study  of  the  cross  sections  for  electron-impact  excitation  of  the  low-lying 
excited  states  of  CO  from  threshold  to  about  4  eV  above  threshold.  These 
studies  should  provide  the  most  reliable  estimates  of  the  cross  sections  for  ex¬ 
citation  of  CO  in  the  technically  important,  but  experimentally  challenging, 
near-threshold  region.  The  availability  of  both  measured  and  calculated  dif¬ 
ferential  and  integral  cross  sections  offer  a  unique  opportunity  for  an  in-depth 
comparison  of  the  measured  and  calculated  cross  sections.  Cross  sections  were 
also  calculated  for  impact  energies  well  beyond  the  threshold  region. 

•  We  have  also  carried  out  a  combined  computational  —  experimental 
study  of  the  cross  section  for  electron-impact  excitation  at  the  low-lying  ex¬ 
cited  states  of  N2.  Differential  and  integral  cross  sections  were  calculated  for 
excitation  of  seven  low-lying  states  (A  3£+,  B  3IIfl,  W  3An,  a  1tt9,  w  1AU, 
a7  1S~,  and  B7  3£~).  Critical  comparisons  of  these  calculated  cross  sections 
with  recent  measurements  by  S.  Trajmar  of  the  Jet  Propulsion  Laboratory 
were  also  carried  out  for  the  A  3£+,  B  3n5,  W  3AU,  and  a  lIig  states.  These 
studies  have  provided  the  first  robust  values  of  the  cross  sections  for  excita¬ 
tion  of  the  metastable  A  3£+  state  of  N2  at  near-threshold  electron  impact 
energies. 

•  We  have  studied  the  cross  sections  for  near-threshold  excitation  of  the 
lowest-lying  excited  state  of  OH  by  low-energy  electrons.  OH  fragments  are 
produced  in  flows  of  interest  via  dissociation  of  water  vapor.  No  measured 
data  are  available  for  these  cross  sections. 

•  Cross  sections  for  electron-impact  excitation  of  H20  to  low-lying  excited 
electronic  states,  which  subsequently  fragment  to  H  and  OH,  were  calculated. 
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